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ABSTRACT 

The successful application of holography to the study of three dimensional 
flow fields due to phase objects has been reported in the literature. The 
present report extends this technique to the study of density fields around 
opaque bodies as would normally be encountered in wind tunnel experiments. 

The density field around a 10 degree half- angle cone at zero and ten degree 
angle of attack has been investigated by means of the finite fringe holographic 
interferometry. The three dimensional density field obtained from the reduc- 
tion of the interferograms was found to agree with that obtained from an 
analytical solution of the governing equations. 

This report was supported by: Airtask No. A310310A/551A/1 R010-03-01 



TABLE OF CONTENTS 



Introduction 
Appendix I 



INTRODUCTION 



This is a report on holographic investigations directed towards the 
determination of complex three-dimensional density fields. The report 
consists of essentially Appendix I which includes the thesis of R. C. 

Jagota on the determination of density fields about a cone at angle of 
attack in a wind tunnel. 

A paper based on this thesis, "Finite Fringe Holographic Inter- 
ferometry Applied to a Right Circular Cone at Angle of Attack, by R. C 
Jagota and D. J. Collins, has been accepted for publication by the Journal 
of Applied Mechanics. 

Further work on the determination of the density in corner flows was 
begun towards the end of the contract year. This work is directed at determining 
the flow in a wing fuselage junction at transonic speeds. 

A preliminary program using a flat plate- fin in the Naval Postgraduate 
School wind tunnel has been started at Monterey. This preliminary program 
is designed to solve some of the holographic problems associated with the 
transonic test that are to be conducted at NSRDC. Further reports on this 
work will be forthcoming in December. 
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ABSTRACT 



The successful application of holography to the study of 
three dimensional flow fields due to phase objects has been 
reported In the literature. The present report extends this 
technique to the study of density fields around opaque bodies 
as would normally be encountered In wind tunnel experiments. 
The density field around a ID degree half-angle cone at zero 
and ten degree angle of attack has been investigated in the 
Naval Postgraduate School supersonic wind tunnel. In these 
experiments, the finite fringe method for the production of 
interferograms has, for the first time, been applied to 
holographic interferometry. The density field obtained from 
the reduction of the interferograms was found to agree with 
that obtained from an analytical solution of the governing 
equations as reported by D. J. Jones in Reference 1. The 
computer program used for the reduction of the interferograms 
has been evaluated for the effect of the presence of the cone 
and the shock wave and has been found to yield stable and 



accurate results. 
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INTRODUCTION 



I . 

For over three decades, the Mach-Zehnder interferometer 
has been widely used and has come to be accepted as a stan- 
dard tool for the measurement of the flow field around 
bodies. Until the advent of the laser, its use required 
extreme stability from vibrations and precision optical 
components. Even with a laser as a light source, the 
resulting interf erograms are purely two-dimensional in 
nature in that they provide fringe information for only one 
direction of view through the wind tunnel. Furthermore, 
due to the limit imposed by existing camera shutter speeds , 
the interferogram obtained does not truly represent the 
instantaneous flow field around the body but rather is an 
average over the interval of shutter opening. 

The advent of holography together with the development 
of the high power Q-switched laser has introduced a new 
dimension in the field of flow measurement. Apart from 
removing the need for precision optical quality components 
in the case of double exposed holograms, the use of holog- 
raphy has, for the first time, enabled a three-dimensional 
view of the flow field to be obtained from a single hologram. 
The high power of the laser and the Q-switch mechanism have 
together enabled film exposure times in the region of 20 
nano-seconds to be attained thus, in effect, "freezing" the 
flow during the hologram production process. 
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Techniques for the application of holography to inter- 
ferometry have been reported in Reference 1 (Heflinger 
et al, 1966) and in Reference 2 (Brooks et al , 1966 ). 
Holographic interferometry has been used to determine the 
three-dimensional asymmetric field produced by free-jet flow 
(Reference 3). In this report, a series of holograms placed 
in a 90 degree sector around the jet were obtained simul- 
taneously from a single exposure. Interferograms were then 
obtained from the holograms at suitable angular orientations 
within a 90 degree field of view. This provided sufficient 
information to determine the entire flow field since the jet 
was planar symmetric. Reduction of the fringe data to den- 
sity information was accomplished by expressing the density 
in a series of orthogonal polynomials and inverting the 
fringe-density integral equation by a computer program using 
a scheme reported in References 4 and 5. 

In the present report, the applicability of holography 
interferometry has been extended to those three-dimensional 
density fields produced by opaque bodies as would be 
encountered in wind tunnel experiments. The experiments 
were conducted in the Naval Postgraduate School supersonic 

blowdown wind tunnel at a Mach number of 2.8 using a ten 

* 

degree half-angle cone at zero and ten degrees angle of 
attack. In this case, the multiple hologram technique was 
not employed. Since the flow fields investigated were steady, 
the procedure adopted was to obtain holograms to cover var- 
ious sectors of the flow field by rotating the cone within 
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the wind tunnel, this being a method that would be fully 
applicable to any conventional wind tunnel setup. 

Whereas in an infinite fringe field, the fringe infor- 
mation is obtained in a discontinuous manner, that is, at 
the points of maxima and minima, the interferogram from a 
finite fringe field can be evaluated at any desired point 
quite readily. This is of particular significance in the 
case of density fields where the variations in density are 
small. In the present case, the total variation in fringe 
number was generally in the region of only one to two 
fringes. The finite fringe method was also selected since 
it is less sensitive to vibrations than is the infinite 
fringe field, a problem encountered in Reference 3- The 
finite fringe field was produced by a translation of the 
diffusing screen in the scene beam a distance of between 
.025" and .030" between the two exposures of the hologram. 

A vertical fringe field, obtained by horizontal translation 
of the diffuser, was chosen in spite of the greater effort 
required to reduce the interf erograms at a particular section 
owing to the greater uniformity and clarity of the fringes 
thus obtained as compared to horizontal fringes. Fringe 
data obtained from the interferograms was reduced to density 
using the same computer program as used in Reference 3- The 
density field obtained experimentally was compared with the 
analytical solution reported in Reference 6, interpolating 
the results therein for the actual mach number of the experi- 
ments. Good correlation was achieved, thus demonstrating 
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the utility of holographic interferometry for normal wind 
tunnel experimentation. 
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II. EXPERIMENTAL APPARATUS 



A. THE WIND TUNNEL 

This investigation was conducted in the Naval Postgraduate 
School blow-down supersonic wind tunnel. The wind tunnel is 
a blowdown-to-atmosphere facility which has a test section 
4 inches by 4 inches in cross-section and 6 inches in length. 
Interchangeable nozzles are available for producing nominal 
test section mach numbers of 4.0, 2.8, and 1.7. The nominal 
run time is 5 minutes at a mach number of 2.8 with a maximum 
stagnation pressure of about 105 pounds per square inch. 

For this series of experiments, the side walls of the wind 
tunnel were replaced by plexiglas walls two inches thick 
having a refractive index of 1.49 so that a complete view of 
the flow in the test section could be obtained. 

B. THE EXPERIMENTAL SETUP FOR PRESSURE MEASUREMENT 

Figure 1 is a schematic representation of the equipment 

arrangement used for pressure measurement. The static 
pressures on 16 out of a total of 30 orifices on the surface 
of the cone were measured by a Scanivalve (Model 48J4-1065 
with a 0-10 psia. transducer) which was calibrated in inches 
of mercury vacuum. The measurements were recorded on a 
Honeywell 2106 twelve channel Viscicorder. The scanivalve 
was capable of scanning 15 ports per second. Photographs of 
the experimental setup are provided in Figures 2 and 3- The 
plexiglas side walls of the wind tunnel with the Scanivalve 
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mounted on top of the test section are shown in Figure 2. 

The power supply, control box and Viscicorder can be seen in 
Figure 3. 

C. THE HOLOGRAPHIC ARRANGEMENT 

A schematic of the holographic arrangement employed is 
shown in Figure 4 and photographs thereof are in Figures 
5(a) and 5(b). The equipment was mounted on a table that 
rested on a portion of the floor that was vibrationally 
isolated from the rest of the building. The monochromatic 
light source used was a Korad K-l pulsed ruby laser operating 
at a wavelength of 6943 Angstroms together with a Pockels 
cell Q-switching device. The resulting effective exposure 
time was about 20 nanoseconds and, with a cavity length of 
the laser of 37 cms . , a coherence length of about 2 cms . 
was obtained. A side band geometry was employed to obtain 
the holograms of the flow field around the cone. A diffuse 
glass in the path of the scene beam was used to produce 
light field holograms, which because of the diffusion of 
the direction of the rays in the scene beam, produced an 
effective field of view in a hologram of approximately 16 
degrees. A continuous wave helium-neon laser was used for 
alignment of the Q-switched laser and the optics. A Lauda 
constant temperature circulator Model N controlled by an 
electronic relay type R-10 was used in conjunction with a 
Culligan de-ionizer to maintain the laser head and output 
etalon at a constant temperature of 28.4 degrees Centigrade. 
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In Figure details of the side band geometry employed 
can be seen. The reference beam was led on top of the wind 
tunnel by means of the two mirrors M and the size of the 
beam could be adjusted by means of lenses U. The positions 
of the diffuser plate and the two grids, one on each of the 
wind tunnel side walls, are also shown therein. In Figure 
5(a), the Korad K-l laser and the smaller alignment laser 
are shown together with the arrangement of the optical 
equipment. The diffuser plate mounted on the precision X-Y 
slide can also be seen. The holographic table and test sec- 
tion was completely enclosed in a wooden box to enable the 
experiments to be conducted in the daytime. This is shown 
in Figure 5(b). Also shown therein are the Korad K-l laser 
power supply and shutter control assemblies, the Culligan 
de-ionizer and the Lauda constant temperature circulator. 

In order to be able to view the cone at angle of attack 
always in a direction perpendicular to its axis when it was 
rotated in the wind tunnel, it was necessary to construct 
separate grids for each rotational angle at which the holo- 
grams were taken. The front grid used for the 82.5 degree 
rotational position of the cone is shown in position in 
Figure 6. The grid lines can be seen to be parallel and 
perpendicular to the cone axis for this particular angular 
position. The orientation of the grids was designed to be 
at the correct angle for the particular rotational position 
of the cone. Since, on rotation, the axis of the cone was 
not in a vertical plane, it was necessary to offset the two 
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grids with respect to each other in such a manner that on 
lining up the grids, the correct viewing position perpendic- 
ular to the cone axis was automatically obtained. The 
method of obtaining the orientations and the relative offset 
of the grids is shown in Figures 7(a) and 7(b). Since the 
cone spindle was bent at a distance of 0.7 inches from the 
base of the cone, the various rotational positions of the 
cone tip are indicated in the end view in Figure 7(a). The 
corresponding positions of the cone axis in the side and 
bottom views are also shown therein. The orientation of 
the grid lines for each rotational position of the cone as 
shown in Figure 7(b) were obtained from the positions of 
the -cone axis in the side view for each particular case. In 
order to be normal to the cone axis, lines perpendicular to 
the orientation of the cone axis in the bottom view were 
transcribed, each of these representing a line of sight. 
Since the distance between the grids was 7 7/8 inches, the 
necessary amount by which the two grids were required to 
be offset for each rotational position of the cone tip were 
obtained as shown in Figure 7(b). 

D. THE WIND TUNNEL MODEL 

The cone model used is shown in Figure 8. The 10 degree 
semi-apex angle cone model was fabricated out of stainless 
steel and had pressure ports at three cross-sections, the 
angular position of the ports being so arranged as to lie on 
a ray from the apex of the cone. Since a fixed mount was 
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employed, provision was made to attach the cone to spindles 
bent to various angles to achieve the desired angle of attack. 
Stainless steel tubing of .036 inches internal diameter was 
used from the static pressure orifices on the surface of the 
cone up to the base of the model, the diameter of this being 
adequate to provide sufficient time response of the Scani- 
valve. Any fifteen out of the thirty ports could be connec- 
ted at a time to the Scanivalve through the hollow cone 
spindle and mount, this restriction being imposed by the 
internal diameter of the spindle. 

Rotation of the model about the horizontal axis of the 
spindle necessitated drilling 4 holes on each side of the 
wind tunnel side walls through which the model spindle could 
be unlocked from the mount and rotated. The plexiglas side 
walls of the test section and the holes therein can be seen 
in Figure 6. The collar on the model spindle containing 
radial holes (also visible in the figure) was used to 
rotate the cone by means of a spike arrangement inserted 
through one of the holes in the tunnel side walls. The 
angle of rotation was established by lining up graduations 
on the spindle sleeve with a graduation on the mount. The 
model and spindle used for the axisymmetric case can be seen 
in Figure 9. The pressure ports at the three stations and 
one of the pressure tubes from the cone passing through the 
hollow spindle are visible in the photograph. After the 
spindle was screwed into the cone, it was locked in position 
by means of the setscrew shown. 
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III. THEORETICAL CONSIDERATIONS 



A. EVALUATION OF THE DENSITY FROM HOLOGRAPHIC INTEPFFROGRAMS 

1 . The Basic Equation of Interferometry 

The interferometry principle is based on the fact 
that two originally coherent rays of light will reinforce or 
annul each other at a screen on which both rays are pro- 
jected depending on their relative phase there. The phase 
difference is directly a function of the difference in the 
lengths of the optical paths traversed by the two rays. If 
the optical path is changed by an amount NX, where X is the 
wavelength of the light and N is an integer, then the order 
of interference changes by an amount N, that is, a shift of 
N fringes occurs. The change in optical path can be related 
to the density through the index of refraction and the 

Gladstone-Dale constant. If the speed of light in a medium 

c 

is represented by -pp- where c q is the speed of light in 
vacuum and n is the index of refraction, the additional time 
required to traverse the test section due to a change in the 
index of refraction from n^ to n^ can be calculated to be 
equal to At = — (n ? - n, ) and thus the change in optical path 
obtained is c Q At , that is, L(n 2 ~n-^). The fringe shift, g, 
which is equal to Ap is thus given by LCn^-n^). 

For a given substance and for a given wavelength of light, 
the index of refraction is a function of density. In the 
case of gases, since the speed of light varies only slightly 



from that in vacuum, it can be represented very accurately 
by the first two terms of the Taylor series: 



n - / + p £ ( 1 ) 

S 

!>' 

where 

p = reference density at CTC, 760 mm. Hg. 
s 

3 = dimensionless constant related to the Gladstone-Dale 

constant by: 

K = £ 

f s 

The value of 3 for air at A = 5893 Angstroms is 0.00n2 Q 2, 
the variation with wavelength being small. The fringe shift 
relation for a fixed difference in index of refraction 
between the two beams becomes: 



? ■/’H (2 

If the density is variable in the test section, the net 
change in optical path is the integrated effect along the 
ray which gives the relationship: 







U 

$ J f (*> >> *c) elf 

O 



( 3 ) 



where 
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( 5 ) 
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ds = Incremental distance along 1 a ray 
z c = a plane of constant z. 

Equation (3) is the basic equation which has to be 
inverted to obtain the unknown density from the known fringe 
shift values obtained from an interferogram. 

2 . The Integral Inversion 

The integral inversion method used was first repor- 
ted by C. D. Maldonado et al in 1965 (1965; Olsen, 1968) and 
subsequently used by R. D. Matulka (Reference 3) to calculate 
the density variation in an asymmetric free jet from the 
fringe numbers obtained from holographic interferograms . 

The procedure involves the representation of the function 
f(x,y,z c ) by a complete set of orthogonal functions, with 
the unknown coefficients being evaluated by use of the 
orthogonality relationship between the set of functions. 

The set of functions employed have a region of orthogonality 
that covers the entire plane and have the property that they 
are "invariant in form" to a rotation of the co-ordinate 
system. The co-ordinate system used for the inversion is 
shown in Figure 10 where X and Y represent fixed laboratory 
co-ordinates and X 1 and Y 1 the co-ordinates in which the 
fringe number function is defined. These co-ordinates, 
therefore, rotate with respect to the fixed co-ordinates X 
and Y as the angle of view through the test section is 
varied . 

In operator form, equation (3) can be represented as: 

H s, *c) = T , *c ) (6) 
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and the object of the procedure is to evaluate f, that is, 
to obtain: 



J- O 



T ‘ 



H (f, X 




( 7 ) 



where T”^- represents the inverse of T. 

This inversion is achieved by utilizing a pair of 
orthogonal polynomials U^ 2k (ax,ay) and H m+2k (ay r ) which are 
related by the transform relationship: 



T 



I U 1K ** 



) -e / - _e / . H («y) (r; 

1 f~ • i j A\ i '7 y -)♦■** t2K ^ 

J r*U"> r ^) 1 1 * 



’Z. L ^ ^ 

£k! C rr,r ^){ J 



where : 



H («»') • /*'T'S h d(«*f][ s'*' 1 H (**')( 

1-Xk m J 



( Q ) 



and are the Hermite polynomials. 

The function f(x,y,z ) is first expanded in a series of 
the polynomials U m+2k as: 



o *li c * M «) u3r^) e 

•v\ s. c *+\ -r XK 



i 2 l '\ 

~0*x rex y / 



(10) 



where e = Jg for m = 0 and 1 for m = 1,2,3, , ^ is an 

m 3 3 3 5 

+m 

arbitrary scale factor, and C m+2k are the unknown expansion 
coefficients. In order to evaluate these coefficients, an 
expression for the function e(S ,y^z c ) is obtained in terms 
of the polynomials by transforming f(x,y,z ) in accordance 
with equation (6). Using the transform relation given by 
equation (8), the following expression is obtained: 
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The expansion coefficients can then be evaluated by 
applying the following orthogonality relationship to 
Equation ( 11 ) : 



r 

f rt ' 

I * 



f r 

•6 cL 5 



/ /-/ 6</J h.,J 7 p') - 

J .*+ 2 K 



*/T 



JJL*[o » <* 



(12) 



that is, by taking the scalar product of equation (11) with 
T 5 i- / , / \ 

.£ ) and using equation (12), the coefficients of 

the series an equation (10) are obtained as: 

<7 W * - J 'titled 7 K) (13) 

l>r 't i 

Substitution of equation ( 13 ) into equation ( 10 ) results 
in the density variation being given by: i 

/^-w 7 zz ^ 



+ 

The functions employed are complex polynomials 



defined as : % r- , 7 > 

- P- >• 4 ■ 
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1 ■ ^ 

where <f> = tan ^-r and L™ are the Associated Laguerre 
polynomials : 
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Substitution of equation (15) into equation (1^1) yields: 
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where : 
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Equations (17), (18), and (19) are the basic equations 
used to obtain the density distribution of a completely 
asymmetric flow field. 

3 . The Numerical Procedure 

The calculation of the unknown coefficients B ir . n , (a) 

m+2k 

and D^ + 2 ^(a) in the series representation of f(x,y,z ) as 
given by equation (17) is not analytically possible in gen- 
eral since the function g(y',£,z c ) representing the variation 
of fringe number is an experimentally measured quantity and 
thus not explicitly given as a function of the co-ordinates 
of a field point . If the analytical procedure is to be 
used to obtain the density from experimental observations of 
the fringe numbers, therefore, it is necessary to evaluate 
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the double integrals in equations ( 1 8 ) and (19) numerically. 
This is achieved by first noting that for every g(y r ,£,z ) 
encountered experimentally, there is a circular domain out- 
side which the fringe distribution is zero so that it is 
possible to replace the limits of integration of +°° and 
in equations ( 1 8 ) and (19) by finite values. The function 
g(y',£,z ) is then approximated by constant step values over 
elemental areas into which the domain is partitioned. This 
results in the B and D coefficients being represented by the 
double series: 



f- 



*■/ ri 
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Using the derivative formula for Hermite polynomials, 
the following equations result: 

B M I t 9(Sj 
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Since it is not possible to sum over an infinite number 
of terms, equation (17) is replaced by the sum of a finite 
series : 







CcS (m cp') 



yj> 

'*»*> t Zh 




The accuracy of the representation depends on the values 
of K, M, A E, and Ax chosen and convergence of the series is 
influenced by the parameter a. 

B. THE FLOW AROUND A RIGHT CIRCULAR CONE 
1 . Axisymmetric Case 

The problem of the supersonic flow around a cone has 
been extensively studied by several investigators ever since 
it was first opened up by Busemann, Bourquard and Taylor and 
Maccoll. The basic approach to the problem is outlined below 
but a detailed description may be found in References 12, 13, 
14 and 1?. 

Starting with the equations of motion of inviscid 
flow in three dimensions expressed in spherical co-ordinates: 

U. v ']£ Ju. 

d/i ^ 
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do 





( 1 ) 
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and the equation of continuity: 



A. (f u h A. (f'l'ir i- A (fsi^ ) £ 

j/l. ^<4 



( 2 ) 



the cylindrical symmetry of the cone allows one to eliminate 
the velocity component w and all terms dependant on <f). In 
addition, application of the condition that the flow be 
conical, that is, no variation of the flow parameters with 
r, reduces equations (1) to: 

cLu. _ ir - O 

"d-G 

Av +. u. + .jL Af - O 
ctC C ctC 

the first of these equations asserting that the flow must be 
irrotational . The equation of continuity similarly reduces 
to : 

d_ ( f ir £■ ) t- £ f d c - C ( 4 ) 

do 

Elimination of v between (3) and (4) yields: 

| -r CL l ofu 

t cl 0 ~ j c l & 

Setting : 

dj2 - dr dr 

Ct.O d f d £ 

the following ordinary differential equation is obtained: 
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From the Bernoulli equation : 



f dt 
J r 






±_ C 



( 8 ) 



and with assumption of a perfect gas satisfying the equation: 

k 



t 

r 



f 

the velocity of sound a can be expressed in terms of the 
velocity components as: 

2 - L ( t-<) ( c* - - v* ) 



(9) 



a 



( 10 ) 



The density, pressure and temperature are then found to 
vary as : 
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For the solution of the governing equation (7) to be 
physically meaningful, it is necessary to satisfy the boun- 
dary condition of no normal velocity at the cone surface 
and the jump conditions at the shock wave. The physical con- 
ditions on either side of a shock wave are given by the Ran- 
kine Hugoniot relations: 
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where, referring to Figure 11, 

a = Angle between the normal to the shock wave and the 
streamlines in front of it. 

B = Angle between the normal to the shock wave and the 
streamlines immediately behind it. 

<p = Angle between the streamline at any point and the 
radius vector from the vertex of the cone. 



U = Freestream velocity. 



a = Velocity of sound in the freestream. 



x 

y 



P /P 
K w ^ 



1 



X 

z 



p /p 
^w ^ 



o 



P-^ = Pressure of the undisturbed gas in front of the 
shock wave. 



p = Pressure immediately behind the shock wave, 
w 

P q = Stagnation pressure behind the shock wave. 



The boundary condition at the shock wave thus reduces 
to the requirement that the ratio (x/z) defined by the 
Rankine Hugoniot relations be equal to the pressure ratio 
(p/P Q ) defined by equation (11a). This together with 
equations (12) defines the physical characteristics of the 
shock wave . 

Many schemes have been employed in order to numer- 
ically integrate the governing differential equation and 
to satisfy the boundary conditions to yield the complete 
solution of the flow field. Results of some of these are 
tabulated in References 6 and 1 4 . For the purposes of this 
report, the more recent results in Reference 6 have been 
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used . 



2. Plow Around Cones with Small Yaw 



The governing differential equation for the flow 
around a yawed cone Is developed on similar lines as that 
for axisymmetric flow with the difference that In this case, 
the velocity components, pressure and density are expressed 
as a series representation about the axial flow values. A 
detailed description of the procedure may be found In 
References 13 s 15 and 17- An examination of the results 
reveals that although the flow is no longer irrotational , 
the flow field can still be regarded as conical. Further- 
more, the shock wave in this case is still a cone of the 
same semi-apex angle as in the no-yaw case with the dif- 
ference that the shock wave is itself yawed at an angle 
which is proportional to the angle of yaw of the cone with 
respect to the freestream. 

In Reference 6, the method employed to obtain the flow 
properties is basically one of iteration. An estimate is 
made of the shape of the shock wave and this shape is 
improved in such a way that the normal velocity at the sur- 
face of the cone is made close to zero. The results are 
presented in the co-ordinate system shown in Figure 12. 

Values of u, v, w, and p at various azimuthal angles <f> 

^00 

are given for various non-dimensional distances between the 
body and the chock wave. Such results are tabulated for 
cones of different half-cone angles, Mach number and angle 
of attack at which the flow can no longer be considered 
conical. For the case of a 10 degree half-angle cone. 
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this limiting angle of attack is 10 degrees at a Mach 
number of 3.0 and 11 degrees at a Mach number of 2.0. 
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IV. DESCRIPTION OF THE EXPERIMENT 



A. NUMERICAL SIMULATION OF THE EXPERIMENT 

The suitability of the computer program HOLOFER for 
handling cases with discontinuities in slope had been inves- 
tigated by means of various test function inputs into the 
program in Reference 3« The accuracy of the inversion in 
these cases was reported to be within 3-8 percent. In order 
to verify the applicability of the program to the present 
investigation, known values of the density along one line 
of sight, obtained from Reference 6 at a particular cross- 
section in the flow, were fed in as input and Mode 1 
operation of the computer program utilized. In this mode, 
the program calculates the fringe number array that would 
result from the specified density distribution and then 
uses this data as the input for inversion to obtain the 
density distribution. One of the primary objectives of this 
investigation was to determine the effect of opaque objects 
introduced into the flow field and the presence of the shock 
wave on the resulting inversion. It was found that the 
effect of the presence of the cone could be reduced by the 
introduction of a fictitious density distribution in the 
area occupied by it of a constant value equal to the density 
at the cone surface. 
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B. EXPERIMENTAL PROCEDURE 



1 . Laboratory Techniques 

It has been shown in Reference 10 that the trans- 
mitted intensity from a diffuse glass falls off below 30 
percent of the incident intensity beyond viewing angles of 
±8 degrees. Thus this represents a limit on the field of 
view obtainable from a single hologram. For the flow around 
a right circular cone at incidence, there is planar symmetry 
and hence a 90 degree field of view is required for complete 
coverage of the total flow field. This coverage can be 
obtained by the simultaneous production of a number of 
holograms in various azimuthal positions as employed in 
Reference 3 or by a rotation of the viewing angle between 
each production of a single hologram. Taking advantage of 
the steadiness of the flow fields investigated, holograms 
covering various sectors of the flow field were obtained 
by rotating the cone to six positions in a 90 degree sector, 
one hologram being obtained in each case. 

Complete angular coverage of the flow field was 
obtained by a rotation of the cone and spindle arrangement 
in the mount by an angle of 15 degrees between each hologram. 
Since it was desired to obtain the flow field at the same 
cross-section for the various angles of rotation, it was 
necessary to obtain the holographic interferograms in such 
a way that the line of sight was always perpendicular to 
the axis of the cone. For this purpose, a pair of grids 
were used, one on each wall of the test section, in which 
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the lines were designed to be parallel and perpendicular 
to the cone axis for a particular rotational position of 
the cone. The two grids were offset laterally with respect 
to each other in order to account for the rotation of the 
cone axis away from the central plane when the cone was 
rotated. On production of the holograms, the variable 
viewing property of light field holograms was utilized to 
line up the two grids thereby ensuring that the line of 
sight was normal to the axis of the cone. This would have 
involved a considerable amount of complexity and effort if 
a conventional Mach-Zehnder interferometer had been employed 
for the same purpose . 

Since the coherence length of the pulsed laser was 
only about 2 cms . at the output used, it was necessary to 
maintain the optical path lengths in the reference and scene 
beams nearly equal. Since the scene beam path traversed 4 
inches of plexiglas while the reference beam passed through 
two thick plano-convex lenses, it was necessary to compen- 
sate for this by making the reference beam path 3-5 inches 
longer than that of the scene beam. A reference to scene 
beam ratio of 4:lvas found to yield good holograms. 

The holograms were obtained on Agfa-Gevaert 8F.-75 
holographic plates, 4" x 4" in size. Development was for 5 
minutes in Kodak D-19 developer, thirty seconds in an acetic 
acid stop bath of standard dilution, five minutes in standard 
fixer, one minute water wash followed by immersion in a 
Kodak PhotoFlo wetting agent prior to drying. Reconstruction 
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of the scene was achieved by illuminating the hologram with 
the light from a helium-neon laser operating at a wavelength 
of 6328 Angstroms. 

2 . Photographic Techniques 

To obtain reconstruction of the image from the holo- 
gram, various methods can be employed. In the usual method 
where the hologram is re-illuminated by a beam similar to 
the original reference beam as shown in Figure 13, each 
point on the photograph is produced by a series of non- 
parallel rays that appear to originate from various sources 
on the diffuse glass that was used to construct the hologram. 
An almost parallel set of rays can be selected by the 
positioning of a small aperture at the focal plane of the 
imaging lens as shown in Figure 14. A similar effect is 
achieved by illuminating the hologram by a conjugate beam 
of small diameter as depicted in Figure 15. Because of the 
small diameter of the reconstruction beam, the illuminated 
portion of the hologram represents a small aperture and thus 
the image produced has a large depth of field. This was of 
considerable advantage in the production of the interfero- 
grams from the holograms since it enabled the front and rear 
grids, the cone and the fringes to be simultaneously pro- 
jected on a screen and thus made the task of obtaining the 
correct viewing angle and the interf erogram quite easy to 
achieve. A photograph of the scene was obtained by focusing 
the camera at the desired plane in the test section, usually 
in such a manner so as to be in the plane of the fringes . 
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The lines of sight recorded on the photograph in this case 
represent the almost parallel pencil of rays from the 
diffuser that pass through the effective aperture in the 
hologram . 

3 . Data Reduction 

To obtain photographic interferograms , the hologram 
was illuminated by the conjugate reference beam using a 
continuous wave laser operating at a wavelength of 6328 
Angstroms. A camera back with viewing screen was placed in 
the position of the real image as shown in Figure 13. By 
lining up the appropriate horizontal and vertical lines of 
the front and rear grids, the viewing angle corresponding 
to the particular line of sight desired was achieved. As 
the spacing of the lines on each grid was ^ inch and the 
two grids were physically 7 7/8 inches apart, a viewing 
angle of 3° 38’ was obtained by lining up a line on one of 
the grids with the adjacent one on the other. A f 7 • 7 lens 
with full aperture and an exposure time of 4 seconds was 
used to obtain workable holograms on Poloroid Type 55 P/N 
film. The reduction of fringe shift was obtained by project- 
ing the negative in a photoenlarger and by tracing the 
fringes, the cone surface, the shock wave and the grid line 
at the desired station at which the reduction was to be 
performed, on a sheet of paper. Some judgement was usually 
necessary in locating the centre of the fringes but in 
general the light fringes were found to be narrower than the 
dark fringes and a location of the fringes to within ±0.1 of 
the fringe spacing in the freestream could be obtained. The 
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fringes in the region outside the shock were extended towards 
the cone surface and the fringe shift read at the points of 
intersection of the displaced fringes with a number of rays 
from the cone vertex. The assumption of locally conical 
flow in the region within two to three fringes on either side 
of the section was made in order to obtain the fringe dis- 
placements at the section from those at the points of 
intersection of the rays from the cone vertex with the 
fringes on either side of the section. The distance between 
the cone tip and the section at which reduction was per- 
formed was measured from the projected image and, since the 
actual distance between these points was known, it was 
possible to calculate the factor by which the projected 
image was magnified as compared with its true size. This 
magnification factor was then used to determine the actual 
radial distances from the cone axis at which the fringe 
shifts were obtained. The radius of the inversion circle 
was selected such that the maximum distance of the shock 
from the cone axis was at 95 percent of the radius of the 
inversion circle. From the data so obtained, the radial 
variation of fringe number was plotted and a smooth curve 
drawn through the points. The fringe number at 201 equidis- 
tant points was then read off from the curve and utilized 
as the input into the computer program HOLOFER in Mode 3. 

This mode utilizes raw data taken at the proper interval 
and directly read into the G array in the program by Sub- 
routine READ. Further details about the use of this computer 
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program are outlined in Appendix C. For the axisymmetric 
case, it was necessary to feed in the fringe data along one 
line of sight only and to obtain the density variation along 
one radial line. For the asymmetric case, fringe data was 
fed in along six lines of sight in a 90 degree field of view 
and the inverted density field obtained along 9 radial lines 
spanning a 180 degree field of view on one side of the axis 
of planar symmetry of the flow. A typical reduction of an 
interf erogram to obtain the radial variation of fringe 
number at a section is shown in Appendix A. 
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V. EXPERIMENTAL RESULTS AMD DISCUSSION 



A. DENSITY FIELD OBTAINED EXPERIMENTALLY 
1 . Computer Simulation Results 

The density distribution obtained from the AGARD 167 
tables for the axisymmetric case was used as a test function 
in Computer program HOLOFER. The fringe distribution cal- 
culated from this density distribution was inverted to obtain 
the original density distribution. A plot of the original 
density distribution and that inverted by the program is 
shown in Figure 16. Since the Mach number for this inversion 
was 2.84, it was necessary to first interpolate the 13 values 
obtained from Reference 6 for a Mach number of 2.84, and 
then to obtain from these, the values at 201 equidistant 
points along a line of sight within the region of the 
inversion circle. This was done using Computer Program 1. 

As a side benefit obtained from a use of this program, the 
discontinuities in slope in the density distribution at the 
cone surface that would otherwise have resulted by assuming 
a constant value in the region occupied by the cone, were 
eliminated. Various values of K, M and interval size were 
experimented with to obtain the minimum deviation of the 
inverted density distribution from the input value in the 
region close to the shock wave . Increasing the value of K 
above 500 was found to have little effect on the accuracy, 
the maximum variation of the inverted density being 1.5 
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percent of the input density value at the point of inversion. 
The value of a was found to affect the total number of terms 
in the series representation of the B and D coefficients 
required for convergence within the value of e specified. A 
value of K = 2000 , M = 5 5 ot = 2.5 and an interval size of 
0.01 was found to yield the minimum effect of the shock 
discontinuity . 

2 . Experimental Results 
a. Axisymmetric Case 

A photograph of the holographic interf erogram 
obtained is shown in Figure 17(a). The fringe data obtained 
from a reduction of the two sides of the cone is shown in 
Figure 17(b). Since the flow was axisymmetric, an average 
of the two readings was taken and a curve plotted. The 
values at 201 equidistant points was read off, this being 
input into Subroutine READ of computer program HOLOFER as 
the G array. The density distribution obtained from the 
output of the inversion is shown in Figure 18. Also shown 
plotted is the density field obtained from the AGARD 167 
tables interpolated for a Mach number of 2.84 by computer 
program 2, this being the Mach number of the tunnel obtained 
from pressure measurements. The pressure data obtained 
from the viscicorder is shown in Figure 19 . The value of 
Mach number calculated from these pressure measurements com- 
pares favorably with the value obtained in earlier experi- 
ments from shadowgraphs of the flow in the tunnel. A 
comparison of the density distributions from the AGARD tables 
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and those obtained from the experiment show good agreement , 
the maximum variation in overall density being 8 percent of 
the theoretical density value, 
b. Asymmetric Case 

The views along which holograms were taken for 
the cone at 10 degrees angle of attack are shown in Figure 
20 and photographs of the interf erograms obtained are shown 
in Figures 21 to 26. Except for the 67.5 degree and 82.5 
degree angles of view, the perturbation to the flow on the 
leeward side of the cone was insufficient to enable either 
the shock position or the fringe displacements on this 
portion of the flow field to be obtained. The fringe shifts 
in this region were, therefore, taken as being zero. A 
comparison of the shock wave position obtained experimental- 
ly from the interferograms with that from the AGARD 167 
tables interpolated for a Mach number of 2.84 is shown in 
Figure 27* The fringe data were input into computer program 
HOLOFER starting with the 82.5 degree fringe distribution 
in order to be consistent with the plane of symmetry for 
which that program was written. The resulting density 
distribution output from the program was along the lines of 
sight indicated in Figure 28 . A comparison of the experi- 
mentally obtained radial distributions with those from 
the AGARD 167 tables interpolated for a Mach number of 2.84 
are shown in Figures 29 to 35. 
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B. ERROR ANALYSIS 



The errors in the final solution are mainly due to errors 
in the fringe data input into the computer program. One of 
the major sources of error was due to the difficulty in 
obtaining the slope of the fringe lines in the freestream to 
a very high degree of accuracy. From the results of the 
experiments conducted on a 10 degree half angle cone in 
Reference 18, it was apparent that the flow in regions close 
to the tip of the cone would be affected by the boundary 
layer on the cone whereas sections near the base of the cone 
would be affected by the wake and expansion of the flow. In 
order to obtain a flow which was locally nearly conical, it 
was thus necessary to reduce the interf erograms at a section 
approximately 6 0 to 70 percent of the cone length downstream 
of the tip. At this section, however, the shock wave was 
about 1 inch from the axis of the cone and considering that 
the fringes in a region about \ inch adjacent to the tunnel 
walls were affected by the boundary layer, this left a region 
of about 3/^ inch in which the fringes of the flow in the 
freestream were present. It was thus difficult to obtain 
the slope of the fringes very accurately and it is estimated 
that errors of approximately ±1 mm. in the fringe data along 
any fringe could have resulted. However, since the fringe 
data along the section was obtained by interpolation between 
different fringes at varying radial distances from the cone 
axis, this would have resulted in a smoothing of the errors 
to about ±0.5 mm. Since the maximum fringe shifts obtained 
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were in the region of 8 mm., an error of ±0.5 mm. corresponds 
to an error of approximately ±8 percent in the magnitudes of 
the fringe shifts obtained. 

The errors between the experimental and theoretical 
density curves in Figure 18 reflect the above estimates of 
error quite well. The experimentally obtained density 
distribution is generally about 5 % higher than the theoreti- 
cal curve and this rises to a maximum of 8 % close to the 
cone surface. Similar errors are apparent in general in 
Figures 29 to 35. However, in these curves the maximum error 
is somewhat higher probably due to the fact that reduction 
was performed at a section closer to the base of the cone 
than for the axisymmetric case. In so far as the general 
magnitude of the density is concerned, the experimentally 
obtained distribution follows the theoretical one fairly 
well. Except for the density distributions for <f> = 0° and 
<J> = 22.5°, the density starts to rise above the freestream 
value at a normalized radial distance from the cone axis 
of about 0.7, then assumes a nearly constant value before 
rising steeply again at a radial distance closer to the 
cone axis . This appears to be inconsistent with the data 
input into the program in which the radial distance of the 
shock from the cone axis is seen to decrease with increasing 
values of <p . 
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V . SUMMARY 



A. CONCLUSIONS 

The finite fringe method for the production of holo- 
graphic interferograms has been applied successfully to 
the determination of the three-dimensional density distri- 
bution of the flow around a 10 degree half angle cone. The 
method has been demonstrated to yield results within an 
accuracy of 8 percent for the axisymmetric case and within 
13 percent for the asymmetric case. The magnitude of the 
errors are mainly a result of the constraints imposed by 
the existing wind tunnel size and the fact that the cone 
model used produced perturbations in the flow field of a 
fairly small magnitude so that small errors in measurement 
were reflected as relatively large percentage errors. The 
computer program H0L0FEP has been evaluated for the effect 
of the shock and the cone and found to yield good results. 
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FIGURE I. SCHEMATIC REPRESENTATION OF THE EQUIPMENT 

ARRANGEMENT USED FOR PRESSURE MEASUREMENT 







FIGURE 2. WIND TUNNEL TEST SECTION AND PRESSURE 
MEASURING EQUIPMENT 
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FIGURE 3. PRESSURE RECORDING EQUIPMENT, CONTROL 
BOX, AND POWER SUPPLIES 
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OF THE HOLOGRAPHIC 
ARRANGEMENT 






FIGURE 5a. HOLOGRAPHIC TABLE LASER CONTROL BOX 
AND TEMPERATURE CONTROL EQUIPMENT 
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DETAILS OF THE HOLOGRAPHIC ARRANGEMENT 
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FIGURE 5b. 




FIGURE 6. VIEW OF TEST SECTION SHOWING DETAILS OK 
MODEL, MOUNT, AND GRID 
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FIGURE 7a. POSITIONS OF THE CONE AXIS FOR VARIOUS ROTATIONAL POSITIONS 
OF THE CONE TIP 




82.5° 67.5° 52.5° 37.5° 22.5° 7.5° 




FIGURE 7b. ORIENTATION OF GRID LINES AND OFFSET OF 
THE GRIDS 
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ALL RADIALLY DRILLED HOLES .040 DIA 
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FIGURE 8. DETAILS OF THE CONE MODEL 





FIGURE 9. THE CONE AND SPINDLE ARRANGEMENT 
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FIGURE 10. CO-ORDINATE SYSTEM USED FOR THE INVERSION 
OF FRINGE NUMBER TO DENSITY 
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Shock Wave 




FIGURE II. CO-ORDINATE SYSTEM USED FOR THE GAS 
DYNAMIC EQUATION 
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SHOCK WAVE 
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FIGURE 12. CO-ORDINATE SYSTEM USED IN REFERENCE 




FIGURE 13. EFFECT OF APERTURE SIZE FOCUS PLANE 

POSITION ON THE PENCIL SIZE OF RAYS ABOUT 
A LINE OF SIGHT RECORDED BY CAMERA 
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FIGURE 14. SPATIAL FILTERING TECHNIQUE FOR SELECTING 
PHOTOGRAPH OF CONSTANT ANGLE LINES OF 
LIGHT 
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FIGURE 15. LENSLESS PHOTOGRAPHIC TECHNIQUE USING A 

CONJUGATE REFERENCE BEAM OF SMALL DIAMETER 
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FIGURE 16 

COMPARISON OF THE RADIAL DENSITY DISTRIBUTION OBTAINED FROM MODE I 
WITH THE AGARD 167 RADIAL DENSITY DISTRIBUTION - AXI-SYMMETRIC CASE 




FIGURE 17a. HOLOGRAPHIC INTERFEROGRAM OBTAINED FOR 
THE AXI- SYMMETRIC CASE 
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Fringe Distance 
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Non-Dimensional Distance from Cone Axis 



FIGURE 17b. RADIAL FRINGE DISTRIBUTION FROM A REDUCTION OF 
THE HOLOGRAPHIC INTERFEROGRAM FOR THE 
AX I- SYMMETRIC CASE 
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FIGURE 18. 

COMPARISON OF THE RADIAL DENSITY DISTRIBUTION OBTAINED FROM 
EXPERIMENT WITH THE AGARD 167 DISTRIBUTION - AXI- SYMMETRIC CASE 
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FIGURE 19. PRESSURE TRACE OBTAINED FROM THE 
FOR THE AXI- SYMMETRIC CASE. 
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FIGURE 20. VIEWS ALONG WHICH HOLOGRAMS WERE OBTAINED 
FOR THE ASYMMETRIC CASE 
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FIGURE 21. HOLOGRAPHIC INTERFEROGRAM FOR 7.5° 
ANGLE OF VIEW 
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FIGURE 22. HOLOGRAPHIC INTERFEROGRAM FOR 22.5° 
ANGLE OF VIEW 
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FIGURE 23. HOLOGRAPHIC INTERFEROGRAM FOR 37.5° 
ANGLE OF VIEW 
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FIGURE 24. HOLOGRAPHIC INTERFEROGRAM FOR 52.5° 
ANGLE OF VIEW 
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FIGURE 25. 



HOLOGRAPHIC INTERFEROGRAM FOR 67.5° 
ANGLE OF VIEW 
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FIGURE 26. HOLOGRAPHIC INTERFEROGRAM FOR 82.5° 
ANGLE OF VIEW 
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FIGURE 27. COMPARISON OF THE POSITION OF THE SHOCK WAVE 
OBTAINED EXPERIMENTALLY WITH THAT FROM AGARD 
167 FOR M = 2.84 , ASYMMETRIC CASE 
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Axis of Planar Symmetry 




FIGURE 28. LINES OF SIGHT ALONG WHICH THE DENSITY 
FIELD WAS OBTAINED AFTER INVERSION, ASYM- 
METRIC CASE 
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AGARD Distribution for M = 2.84 




FIGURE 29. COMPARISON OF THE EXPERIMENTAL RADIAL DENSITY DISTRIBUTION 
WITH THE AGARD 167 DISTRIBUTION FOR THE ASYMMETRIC CASE 
AT M =2.84, <f> = 0° 
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FIGURE 30. COMPARISON OF THE EXPERIMENTAL RADIAL DENSITY DISTRIBUTION 
WITH THE AGARD 167 DISTRIBUTION FOR THE ASYMMETRIC CASE 
AT M= 2.84. < 1 > = 22.5° 
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FIGURE 31. COMPARISON OF THE EXPERIMENTAL RADIAL DENSITY DISTRIBUTION 
WITH THE AGARD 167 DISTRIBUTION FOR THE ASYMMETRIC CASE 
AT M=2.84. <f> = 45° , 



||.6 % Error ° AGARD 167 Distribution for M=2.84 

+ Experimental Distribution 




FIGURE 32. COMPARISON OF THE EXPERIMENTAL RADIAL DENSITY DISTRIBUTION 
WITH THE AGARD 167 DISTRIBUTION FOR THE ASYMMETRIC CASE 
AT M = 2.84, <p = 67.5° 
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FIGURE 33. COMPARISON OF THE EXPERIMENTAL RADIAL DENSITY DISTRIBUTION 
WITH THE AGARD 167 DISTRIBUTION FOR THE ASYMMETRIC CASE 
AT M = 2.84, d> = 90° 



AGARD 167 Distribution for M=2.84 
Experimental Distribution 




76 



FIGURE 34. COMPARISON OF THE EXPERIMENTAL RADIAL DENSITY DISTRIBUTION 
WITH THE AGARD 167 DISTRIBUTION FOR THE ASYMMETRIC CASE 
AT M = 2.84, d> =112.5° 



o AGARD 167 Distribution for M=2.84 
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FIGURE 35. COMPARISON OF THE EXPERIMENTAL RADIAL DENSITY DISTRIBUTION 
WITH THE AGARD 167 DISTRIBUTION FOR THE ASYMMETRIC CASE 
AT M = 2.84, d> = 135° 



Avg. Fringe Dist. = — =2.775 
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FIGURE 36. TRACING OF THE PROJECTED IMAGE OF THE HOLOGRAPHIC 
INTERFEROGRAM FOR THE AXI-SYMMETRIC CASE 
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FIGURE 37. PLOT OF A TYPICAL FRINGE SHIFT VARIATION ALONG A RAY TO 
OBTAIN THE FRINGE SHIFT AT THE SECTION 



TABLE I 



VARIATION OP FRINGE SHIFT WI m H DISTANCE 
ALONG RAYS FROM THE CONE VERTEX 

x = Fringe shift (mm.) 

y = Distance from cone vertex along a ray (ins.) 
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TABLE II 

RADIAL DISTANCE FROM THE CONE AXIS OF THE 
INTERSECTION OF RAYS WI'T’H THE SECTION 



Distance from Normalized distance 

cone axis from cone axis 
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TABLE III 



CALCULATION OF THE RADIAL VARIATION 
OF FRINGE NUMBER AT THF SECTION 



fringe spacing in the freestream = 2.778 mm. 



Normalized distance Fringe shift 
from cone axis (mm.) 
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APPENDIX A 



REDUCTION OF AM INTFRFEFOGRAf TO OPTAIN FRINGE SHIFT DATA 

The reduction of only one side of the interferogram ob- 
tained for the axisymmetric case is indicated as an illus- 
tration of the procedure employed. After projection of the 
negative, the cone surface, the shock waves, the grid line 
at the section concerned and the fringes from the projected 
image were traced on a sheet of paper as shown in Figure 3^. 
The distance between the cone surface and the shock at the 
section was divided into a convenient number of parts and 
rays drawn from the tip of the cone through these points. 

The fringe shift distances at the point of intersection of 
each displaced fringe and the rays were measured by means of 
a 7X PEAK scale magnifier as well as the average fringe 
spacing in the freestream. m he radial distance from the 
cone axis to the intersection of each of the rays with the 
section was also measured as also the distance along each 
ray to the intersection with each fringe and with the 
section concerned. A tabulation of these results is shown 
in Tables 1, 2, and 3. For each ray, the fringe shift dis- 
tance was then plotted against the distance along the ray 
as in Figure 37 , and from these curves the fringe shift 
distance at the section was obtained. These fringe shift 
distances were divided by the average fringe spacing in the 
freestream to obtain the fringe numbers at the various 
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points concerned. From a measurement of the distance from 
the tip of the cone to the section on the projected image 
and comparison with the known distance, the magnification 
of the image with respect to the actual conditions existing 
in the wind tunnel was determined. The radius of the inver- 
sion circle was obtained by assuming the distance of the 
shock from the cone axis to be at 95$ of the radius of the 
inversion circle. The distance from the cone axis of the 
various intersections of the rays from the cone vertex and 
the section were obtained as a fraction of the inversion 
circle radius. The fringe number at these points was then 
plotted against the non-dimensional distance from the cone 
axis as shown in Figure 17 and a smooth curve drawn through 
these points. The value of the fringe number in the region 
occupied by the cone was taken as a constant value equal to 
that at the cone surface. From the curve so obtained, 
values of the fringe number at 101 equidistant points were 
obtained on each side of the cone axis (including the value 
at the axis) so that a total of 201 data values resulted. 
This fringe data was read into Subroutine RFAD of computer 
program HOLOFER and inverted to obtain the required density 
field. 

The use of 201 data points was essentially dictated by 
the necessity to be able to define the shock position 
accurately and to have the fringe distribution in the region 
between the shock and the cone described fairly well. Since 
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the fringe curve from which these points were obtained was 
plotted using only about 10 experimental fringe values, the 
use of 201 points does not imply a higher accuracy in the 
density distribution output from the computer program. 



85 



APPENDIX B 



USE OF COMPUTER PROGRAM "HQLOFEF" 

This computer program is designed to invert the array of 
fringe numbers to obtain the associated density field using 
the inversion first proposed by C. D. Maldonado and described 
in Section III of this report. The computer program can be 
run in basically 3 modes as described below: 

( a ) Mode 1 

This mode provides the basic testing capability of 
the program and uses various test functions listed in 
Subroutine FUNCT in order to generate a G array which is 
then inverted back to obtain the original input function. 
Functions other than those specified in FUNCT can also be 
read in on cards by specifying a test function number of 
8.0 in the original list of 42 parameters read in into the 
main program HOLOVERT. In this case the main program first 
calls Subroutine FREAD to read in the data cards. The first 
card in this data deck consists of two values indicating the 
total number of cards to be read in and the ’ Z ' section at 
which inversion is being performed and is input according to 
format 89 of the subroutine. Thence follows one point per 
card according to format 88 representing the numerical test 
function being input. This mode was employed in the present 
case for inverting the 201 points that were obtained from 
the AGAFD tables by interpolation using Computer Program 1. 
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